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Abstract 

We consider moduli fields interacting with thermalized relativistic matter. We deter- 
mine the temperature dependence of their damping rate and find that it is dominated 
by thermal effects in the high temperature regime, i.e. for temperatures larger than 
their mass. For a simple scalar model the damping rate is expressed through the known 
matter bulk viscosity. The high temperature damping rate is always smaller than the 
Hubble rate, so that thermal effects are not sufficient for solving the cosmological 
moduli problem. 
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Many theories beyond the Standard Model of particle physics, such as string and 
supersymmetric theories, contain directions in field space that have a flat potential in 
the supersymmetric limit and couple to other fields only through Planck scale sup- 
pressed interactions. If these directions are stabilized by the same mechanism which is 
responsible for supersymmetry breaking, a particle whose mass is of the order of a TeV 
and a dangerously large lifetime may result. Such fields are often called moduli. If one 
of the moduli starts at a value displaced from the minimum of the effective potential, it 
will perform coherent oscillations around its minimum. The energy of these oscillations 
decreases like the energy density in non-relativistic matter with the expansion of the 
Universe, p oc a~ 3 , while for relativistic matter p oc a -4 . Therefore the moduli may 
soon contribute significantly to the total energy density of the Universe. 

In the vacuum the decay rate of a moduli field with mass m is of order m 3 /M 2 , where 
M = 2 x 10 18 GeV is the reduced Planck mass pQ. The corresponding long lifetime 
r = ~ 10 5 (TeV/m) 3 s causes the cosmological moduli problem Unless the 

initial amplitude of the moduli is very small, they dominate the energy density of the 
Universe at some point, spoiling the success of standard big bang nucleosynthesis. 

Many solutions to this problem have been suggested. Recently it was pointed out in 
Ref. [7j that in the hot early Universe the moduli oscillations could be damped much 
more strongly than in the vacuum. It was found that the damping is in fact so strong 
that the moduli would adiabatically follow the minimum of the effective potential, 
making their relic density harmless for cosmology In this note we reconsider the 
thermal damping of scalar moduli fields. Unfortunately, we cannot confirm the result 
of Ref. 

A modulus field starts oscillating once the Hubble rate H drops below its mass m, 
which happens at the temperature T osc ~ \JmM ~ 10 10 GeV-^/m/TeV. The thermal 
corrections to the mass have a magnitude <T 2 /M, which is of the same order as 
H if the Universe is radiation-dominated. Thus, unless they come with a very large 
coefficient 8J, they do not affect the parametric size of T osc , and they can be neglected 
for T < T osc . 

Thermal damping effects can become important if the damping rate becomes larger 
than the Hubble rate. In order to see whether this happens, we can neglect the ex- 
pansion of the Universe. We consider temperatures T 3> m, because, as we shall see, 
only in this case does thermal damping dominate over the vacuum damping. Then the 
matter fields, i.e. fields that have unsuppressed interactions, evolve on time scales of 
order T _1 , while the modulus field evolves much more slowly. 
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To describe the time evolution of this non-equilibrium system, we imagine starting 
from thermal equilibrium at t = — oo. For simplicity the zero of the modulus field p 
is chosen such that its expectation value at the considered temperature vanishes. To 
drive the system out of equilibrium we add the term £ ext = jip to the Lagrangian, 
where j is an external c-number field that vanishes for t — > — oo. Then, at some finite 
time, j is switched off again. From then on we have an isolated system which is not 
in thermal equilibrium. The way the system got out of equilibrium does not matter 
because p evolves much more slowly than all other degrees of freedom, since m <C T. 
On time scales over which p changes, the other fields have long come into thermal 
equilibrium and do not remember anything about their past. 

As long as (if) <C M one can use the linear approximation to compute the expectation 
value of p 0: 

(cpix)} ~ J d 4 xA iet {x-x')j(x') (1) 

with the retarded propagator 

A ret (x - x') = (p(x)p(x')) rct =iQ(t- ([<p(x), cpix')})^ . (2) 

Here the expectation value is taken in a thermal ensemble with the initial temperature 3 . 
The Fourier transform of the retarded propagator has an analytic continuation 

into the upper half of the complex frequency plane such that 

A ret (p) = A(p° + ie,p). (4) 

We are interested in spatially homogeneous oscillations and can therefore consider an 
cE-independent j. Without any interaction cp would simply oscillate with the frequency 
m. With interaction, m 2 receives finite temperature corrections given by the real part 
of II. The imaginary part of II leads to exponential damping of the p> oscillations. We 
assume that the damping rate is much smaller than the oscillation frequency, which 
will be verified a posteriori 4 . Then the damping rate is given by 

7 = ImII(m T + ie, 0), (5) 

2rriT 



3 The system may get heated by applying the external field, but this effect does not contribute in 
the linear approximation. 

4 Yokoyama 7 finds a damping rate which is much larger than the frequency. 
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where tut is the thermally corrected frequency. Since tut is much smaller than all 
other relevant energy scales in the problem, we can write 

7 ~ — lim — Imll(u; + ie, 0). (6) 

This limit exists since II is real at zero frequency. Therefore, at leading order in m/T, 
the damping rate does not depend on m 5 . For large T one can neglect the masses of 
the matter particles so that, by dimensional analysis, 

T 3 

7 = f W , (7) 

where / is a function of the (dimensionless) gauge, Yukawa, and scalar couplings. For 
T < T osc we have 7 < mf^/m/M, which is indeed much smaller than m. 

In general the calculation of / is non-trivial. We consider a simple model [7] where 
if is coupled to a massless scalar field Xi which mimics a Standard Model field. We will 
see that for this case the damping rate can be expressed through the bulk viscosity, 
which has been computed elsewhere. The Lagrangian is 

with 



C v = \ (^) 2 - (9) 



= \ (dx) 2 - ^X 4 (10) 

£<PX = £fV(dxf- (11) 

The coupling constant A is supposed to be reasonably large, so that the field x is i n 
thermal equilibrium, but small enough so that perturbation theory can be applied. At 
leading order in M~ l the ip self-energy can be written as 

Tl(p° + ie,p) = / dW- ((d X ) 2 (x) (d X ) 2 (0)> ret , (12) 



5 In Ref. [Jj an m-dependent damping rate was obtained. One cause for this discrepancy appears 
to be that in going from Eq. (21) to Eq. (22) of Ref. |Zj the x-damping rate T p has been neglected in 
favour of m, even though T p 3> m. 



3 



where the 2-point function in Eq. (|12|) has to be evaluated for M — > oo. In this limit 
the operator (dx) 2 is related to the trace of the energy momentum tensor of \ 

T „ v = d»xd v X - V^£ x ( 13 ) 
through the trace anomaly [TT] . 

T^ = -^-(d X ) 2 + ---. (14) 

Here /3(A) = 3A 2 /(167r 2 ) + 0(A 3 ) is the (3 function of the scalar x 4 theory. The ellipsis in 
Eq. (jnj) denotes total derivatives and terms which vanish by the equations of motion, 
neither of which contribute to the imaginary part of IT in the limit (JBJ). Combining 
Eqs. (JHJ), (fT2|). and (|14j) . we obtain at leading order in A 

The right-hand side is related to the bulk viscosity ( of a gas of x particles through 
the Kubo relation ^2] 

C = £ lim -Im / d 4 xe^°(T%(x)T%(0)) ret . (16) 

We thus find that the damping rate of the moduli field <p in this model is directly 
proportional to the bulk viscosity of the plasma 



7 = -^-^C (17) 
7 8M 2 A 2 ^ { ! 

In order to compute £ even at leading order in A one has to sum an infinite set of 
diagrams This turns out to be equivalent to solving an appropriate Boltzmann 
equation, which accounts for particle-number-changing processes. The result is ^3] 

C = , b x, Ahi 2 (£A)T 3 , (18) 
S 6(32vr) 4 ^ J ' y J 

where b = 5.5 x 10 4 and f = exp (15C(3)/vr 2 ) /96 = 0.064736. Then we finally obtain 

4 ln 2 (£A) T 3 , s 

7-4V-^ (19) 

for the damping rate in the model (jSJ). 
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We have thus seen that the damping rate of moduli fields in a thermal environ- 
ment is parametrically larger than in the vacuum, when T > m. However, it is always 
small compared to the Hubble rate. Therefore thermal effects alone cannot solve the 
cosmological moduli problem. 
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